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ABSTRACT 

We present results of a recent Chandra X-ray Observatory observation of the central compact object 
(CCO) in the supernova remnant Cassiopeia A. This observation was carried out in an instrumental 
configuration that combines a high spatial resolution with a minimum spectral distortion, and it 
allowed us to search for pulsations with periods longer than « 0.68 s. We found no evidence of 
extended emission associated with the CCO, nor statistically significant pulsations (the 3a upper 
limit on pulsed fraction is about 16%). The fits of the CCO spectrum with the power-law model yield 
a large photon index, F « 5, and a hydrogen column density larger than that obtained from the SNR 
spectra. The fits with the blackbody model are statistically unacceptable. Better fits are provided by 
hydrogen or helium neutron star atmosphere models, with the best-fit effective temperature kT^ « 
0.2 keV, but they require a small star's radius, R = 4-5.5 km, and a low mass, M < O.SMq. 
A neutron star cannot have so small radius and mass, but the observed emission might emerge 
from an atmosphere of a strange quark star. More likely, the CCO could be a neutron star with 
a nonuniform surface temperature and a low surface magnetic field (the so-called anti-magnetar) , 
similar to three other CCOs for which upper limits on period derivative have been established. The 
bolometric luminosity, LJ^j ~ 6 x 10'^'^ erg s~^, estimated from the fits with the hydrogen atmosphere 
models, is consistent with the standard neutron star cooling for the CCO age of 330 yr. The origin of 
the surface temperature nonuniformity remains to be understood; it might be caused by anisotropic 
heat conduction in the neutron star crust with very strong toroidal magnetic fields. 
Subject headings: stars: neutron — supcrnovae: individual (Cassiopeia A) 



1. INTRODUCTION 

Cassiopeia A (Cas A) is the famous remnant of a 
Type II supernova explosion that was apparently de- 
tected about 330 years ago (Aschworth 1980). The com- 
pact X-ray source at the center of Cas A was discovered 
in the first-light Chandra observation (|Tananbaumlll999f) 
and then found in archival ROSAT and E instein images 
■ (jAschenbachI 119991 : iPavlov fc Zavhnlll999( l. It is consid- 
ered as a prototype of the so-called compact central ob- 
jects (CCOs) in supernova remnants, which are possibly 
neutron stars (NSs) w ith proper ties very different from 

those of radio pulsars (iPaylov et al. 2002). 

The analysis by iPavlov et al.i (|2000f ) and 
IChakrabartv et all (|2001h of first Chandra observa- 
tions of the Cas A CCO with the ACIS and HRC 
detectors has shown that it is a pointlike source with a 
flux Fx « (8-9) X 10"" ergs s'^ cm'^ in the 0.6-6 keV 
band. The source showed no periodicity, with an upper 
limit on the pulsed fraction of 35% for P > 20 ms. Its 
spectrum was found to be consistent with an absorbed 
blackbody (BB) model [A^h,22 = A^n/lO^^cm-^ « 0.8, 
kT = 0.5-0.6 keV, R = (0.3-0. 5)L»3.4 km, where 
DsA = D/3.4kpc is the distance scaled t o the rem- 
nant's estimated distance; iReed et al] Il995j | or with a 
steep power-law (PL) model (A'h 22 ~ 2, photon index 
F = 3-4). 

Based on the X-ray properties of the Cas A CCO, 
IPavlov et al.l and IChakrabartv et all suggested that its 



thermal-like radiation might be interpreted as emitted 
from hot spots at the NS surface, similar to those seen 
in some young radio pulsars and accreting X-ray pul- 
sars. However, this source shows no indications of pulsar 
activity (such as radio or 7-radiation, or a pulsar wind 
nebula), the emitting region is unusually hot compared 
to radio pulsars, and the source flux does n ot show vari- 
ations exp ected for an accreting o bject. IPavlov et al.l 
(|2000( ) and lChakrabartv et all (|2001h speculated that the 
unusual properties of this source might be caused by an 
extremely high magnetic field, B ~ lO'^^'-lO^^ C, similar 
to those of anomalous X-ray pulsars (AXPs) and soft 7- 
ray repeaters (SGRs), commonly known as "magnetars", 
but they concluded that more observations were required 
to understand the true nature of this CCO from the spec- 
tral and timing analysis. 

Since the discovery of the CCO, Cas A has been ob- 
served many times with Chandra. In addition to about 
90 short (1-10 ks) ACIS and HRC cahbration observa- 
tions, two 50 ks ACIS-S3 observations (2000 January and 
2002 February), three 50 ks HRC observations (1999 De- 
cember, 2000 October, and 2001 September), and a 70 
ks ACIS/HETG observation (2001 May) have been car- 
ried out. Finally, Cas A was observed with ACIS-S3 
in 2004 February - May (9 pointings, 40-170 ks expo- 
sures) for a total exposure of about 1 Ms. Unfortunately, 
the vast amount of data collected in those observations 
has not added much to our understanding of the CCO. 
The resuhs of the 50 ks ACIS-S3 observation of 2000 



2 



Pavlov and Luna 



January and 50 ks HRC- S observation o f 200 Octo- 
ber have been reported by iMurrav et al.1 (|2002l ) . They 
found a low-significance period of 12 ms in the HRC-S 
data, which was not co nfirmed in the HRC-S observa- 
tion of 2001 September (jRans om l l2002( ). and concluded 
th at the spectrum wa s con sistent with that r e porte d 
bv iPavlov et all (|2000D and IChakrabartv et all pOOlh . 
iHwang et al.l ( 2004D presented beautiful images and spec- 
tra of the Cas A SNR obtained in the 1 Ms observations, 
but, regarding the CCO, they only mentioned that the 
best single-component model is a BB with kT — 422 ± 6 
eV and R = (0.83 ± 0.03)i:»3.4 km; that fit, however, 
was formally unacceptable (xl = 1-57 for 315 degrees 
of freedom [dof]). The analysis of the ACIS observa- 
tions of 1999-2004 did not show statistically significant 
changes of the CCO flux (,Teter et al..,2004 l. virtually rul- 
ing out the possibility that the CCO emission is due to 
accretion. The spectral analysis of those data has shown 
that the fits with one-component thermal models (BB 
or NS atmosphere) leave large residuals at high energies, 
E > A keV, suggesting that either the surface tempera- 
ture is nonuniform or the emission at higher energies is 
of a nonthermal (e.g., magnetospheric) origin. 

However, the spectral fits of those data, and the source 
properties inferred from those fits, suffered from large 
systematic uncertainties because all the "bare-ACIS" 
(no-grating) observations were taken in full frame mode 
(frame time 3.24 s), which resulted in a significant photon 
pileup^, with ~ 20% pileup fraction for plausible spec- 
tral models. Although such a pileup fraction may look 
moderate, it substantially distorts the spectrum. In par- 
ticular, it leads to an artificial excess of counts at higher 
energies (because two photons arriving within a frame 
are counted as one photon of higher energy), which can 
be easily confused with, e.g., a PL tail expected for active 
pulsars and magnetars. 

Additional systematic errors in the bare-ACIS data of 
1999-2004 were caused by the choice of Graded teleme- 
try format^ (to avoid telemetry saturation from the very 
bright SNR), which made it impossible to correct the 
spectra for Charge Transfer Inefficiency (CTI), and an 
off-axis placement of CCO (to image the SNR around 
the chip center), which blurred the CCO image and ham- 
pered the search for a possible pulsar wind nebula). Fi- 
nally, in all but one of the 1 Ms observations the CCO 
was placed at the boundary between two CCD nodes, 
where "bad pixels" contaminate the data. 

In the ACIS/HETG observation of the CCO, which 
did not suffer from pileup and was telemetered in Faint 
mode, the dispersed CCO spectrum was very strongly 
contaminated by the bright SNR background, while the 
zeroth order image had too few counts and was contam- 
inated by the dispersed SNR image. In addition, there 
were virtually no source counts above 4 keV, where the 
ACIS/HETG effective area is very small. 

The CCO was also observed with the EPIC 
instrument onboard the XMM-Newton observatory 

^ If two or more photons interact with the same detec- 
tion cell within the frame time, they are registered as a sin- 
gle event. This effect is known as pileup (see §6.15 of 
T he Chan dra Proposers' Observatory Guide [POG], ver. 11, 
at |http : //asc .harvard. edu/proposer /POG for a detailed descrip- 
tionj! 

2 See §6.4.12 in POG. 



(|Mereghetti et al.ll2002[ ). However, because of the rel- 
atively low angular resolution, the CCO data were 
strongly contaminated by bright SNR filaments in the 
CCO vicinity. 

Thus, despite the very deep observations taken, sys- 
tematic errors, particularly those caused by pileup, have 
precluded an accurate spectral analysis of the CCO, cru- 
cial for understanding its nature. Therefore, we carried 
out another Chandra ACIS observation of the CCO, us- 
ing a subarray of the only activated ACIS-S3 chip. The 
much shorter frame time of 0.34 s in this configuration 
allowed us to obtain a CCO spectrum virtually undis- 
torted by pileup, and to search for pulsations with a pe- 
riod > 0.68 s. Placing the target very close to the optical 
axis, we were able to image the immediate vicinity of the 
CCO with a high spatial resolution and search for com- 
pact extended emission (e.g., a pulsar wind nebula) that 
might be associated with the CCO. 

In this paper we present the results of this observation. 
The data reduction is described in §2. In §3 we present 
the observational results, including the image analysis 
(§3.1), spectral analysis (§3.2), and timing (§3.3). Impli- 
cations of the results and the possible nature of the CCO 
are discussed in §4. 

2. OBSERVATIONS AND DATA REDUCTION 

We observed the Cas A CCO on 2006 October 19 (Ob- 
sID 6690) with the Chandra ACIS detector for 70.181 ks 
in Timed Exposure mode, using Faint telemetry format. 
The target was imaged on the ACIS-S3 chip close to the 
optical axis (with a standard 20" Y-offset to move the 
target from the node boundary). To minimize the pho- 
ton pileup and improve time resolution, we used a 100 
pixel subarray (8.3' x 0.82' field of view [FOV]) near the 
chip readout and turned off the other ACIS chips. In this 
observational setup, the frame time is 0.34104 s, which 
consists of 0.3 s exposure time and readout (dead) time 
of 41.04 ms required to transfer charge from the image 
region to the frame store region. Therefore, the effective 
target exposure time was 61.735 ks. There were no pe- 
riods of substantially enhanced background during the 
observation. 

We reduced the data using Chandra Interactive Analy- 
sis of Observations (CIAO) software^, ver. 4.1.2 (CALDB 
ver. 4.1.2). For the timing analysis, we transformed the 
event times of arrival to the solar system barycenter us- 
ing the barycen tool. In our image analysis we used 
MARX'' and Chandra Ray Tracer (ChaRT) software^ 
We used XSPEC (ver. 12.0.4) for the spectral analysis. 

3. DATA ANALYSIS AND RESULTS 
3.1. Image 

The 51" X 39" image of the CCO field is shown in Figure 
1. The CCO looks like a pointlike source embedded in the 
SNR background. Using the CIAO wavdetect tool, we 
found the centroid of the CCO image at the coordinates 
a = 23'^23'"27.952" and 6 = -h58°48'42.57" (J2000). The 



^ See http : //cxc. harvard. edu/ciao/ " 

** MARX (Model of AXAF Response to X-rays) is a 
suite of programs designed to enable the user to simu- 
late the on-orbit performance of the Chandra satellite. See 

jhttp : //space . mit . edu/ASC/MARX/] 

The software is available at |http : //cxc ■ harvard ■ edu/ chart/ 1 
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Fig. 1.— Chandra ACIS image, 51" x39", of the Cas A CCO and 
its surroundings in the 0.6—6 keV ba nd, af ter applying the subpixel 
event reposition tool by Mori ct al. {'2001'). The pixel size is 0.123" 
(i.e., 1/4 of the original ACIS pixel size). The inner circle around 
the CCO shows the source extraction region of 1.476" radius used 
for the spectral analysis, while the annulus between 2.46" and 3.94" 
(white circles) is the background extraction region. The CCO is 
surrounded by non-uniform SNR emission, with the average surface 
brightness of 89 counts arcsec"^ in the 2.46"-3.94" annulus. 
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Fig. 2. — Radial profiles of the observed data (black) and the 
MARX-simulated PSF (blue) in the 0.6-6 keV band. The surface 
brigthness was measured in the images with pixel size of 0.0615" 
(1/8 of the original ACIS pixel) in 35 circular annuli with 0.123" 
widths. The simulated PSF profile is for the MARX DitherBlur 
parameter of 0.30". The red histogram shows the sum of the 
sumulated PSF and the model constant background, 89.6 counts 
arcsec"^. 

CIAO celldetect tool yielded a = 23'>23'"27.956" and 
S = +58°48'42.58". The differences of 0.03" and 0.01" in 
right ascension and declination, respectively, exceed the 
formal centroiding uncertainties (e.g., 0.006" per coordi- 
nate for the wavdetect measurement). They, however, 
are much smaller than the uncertainty of the Chandra 
absolute astrometry, « 0.2" for each of the coordin ates 
at the confidence level of 68% (jPavlov et al.ll2009a[ l. as 
estimated from the empirical distribution of radial off- 
sets of X-ray positions with respect to the accurately 
known celestial locations for a sample of point sources^. 
The measured CCO coordinates are consistent with the 
most accurate of the previously measured coordinates 

See § 5.4 and Fig. 5.4 in POG. 



(|Fesen et al.ll2006[) . within the uncertainties. 

To look for an extended component in the CCO image, 
such as a pulsar wind nebula (PWN), and to choose an 
optimal extraction aperture for the spectral and timing 
analysis, we simulated a point source observation using 
ChaRT and MARX and compared the results with the 
observed count distribution. To reach a subpixel spa- 
tial resol ution, we applied the subpixel event reposition 
method (jMori et al.ll2001f ) to both the observed and sim- 
ulated images. For the simulation, we chose the spectral 
model wabs x nsa with the best-fit parameters (see i)3.2l 
and Table [T]) and simulated an ACIS-S3 observation for 
the same position on the detector and the same exposure 
time, 61.735 ks, as those of the actual observation. The 
width of the simulated point spread function (PSF) de- 
pends on the value of the MARX parameter DitherBlur. 
This parameter accounts for the ACIS pixelization and 
aspect reconstruction errors, which may be different for 
different observations. We simulated the PSF for a num- 
ber of DitherBlur values, from 0.20" to 0.40", and found 
that the best match of the simulated PSF to the core of 
the observed image is provided by DitherBlur ~ 0.30". 
Figure [2] shows the simulated PSF radial profile and the 
radial distribution of the detected events up to 4.3" from 
the source centroid. To calculate and plot these profiles, 
we rebinned the observed and simulated images to 1/8 
of the original ACIS pixel and measured the numbers 
of counts in circular annuli with widths of 0.123" (1/4 
of the original pixel size). We see from this figure that 
the observed radial profile exceeds the simulated PSF at 
r > 1", remaining approximately constant at r > 2", 
with the average surface brightness of 89.6 ± 1.5 counts 
arcsec"^ in the 2" < r < 4" annulus. The morphology 
of this extended emission on larger scales (see Fig. 1) 
suggests that it originates from the Cas A SNR (possi- 
bly belongs to a faint SNR filament), i.e., the extended 
emission is not a nebula generated by the CCO. This con- 
clusion is supported by the presence of emission lines in 
the spectrum of this extended emission (see § 3.2), which 
are not expected in the synchrotron spectrum of a PWN. 

The sum of the simulated PSF and the uniform surface 
brightness of 89.6 counts arcsec"^, shown by the red his- 
togram in Fig. [21 is generally very close to the observed 
radial profile within the r < 4.3" circle. It lies slightly be- 
low the observed data (i.e., the observed profile is slightly 
broader than the simulated one) at 0.7" < r < 1.2", but 
the statistical significance of this difference is marginal 
(e.g., the largest discrepancy, in the 7th annulus, is sig- 
nificant at the 2a level). The difference might be caused 
by a nonuniformity of the extended emission component 
at such radii or a minor inaccuracy of the MARX sim- 
ulation, but it is hard to believe that the data excess is 
due to a PWN-like emission in this narrow region. An 
additional support for the pointlike structure of the CCO 
image is provided by the image deconvolution with the 
aid of the arestore script in CIAO (b ased on the Lucy- 
Richardson algorithm), performed by 'Kargaltsev et al] 
12009 ) for the comparison with PSR J1617-5055 that 
is embedded in a PWN. These authors conclude that 
the deconvolved image of the Cas A CCO "preserves the 
pointlike appearance with no extended structure" (see 
Fig. 4 of that paper). Thus, we conclude that our ob- 
servation does not provide any substantial evidence for 
a PWN around the Cas A CCO. 
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Fig. 3. — Signal-to-noise ratio as a function of circular aperture 
radius. 

3.2. Spectral analysis 

As the ceo is embedded in the significant SNR back- 
ground, the source extraction region should be smaller 
than that for an isolated point source. To find an op- 
timal extraction radius, we modeled the dependence of 
the signal-to-noise ratio (S/N) on the radius r of the 
source extraction aperture in the 0.6-6 keV band (Fig- 
ure 3), using the above-described MARX simulation and 
the measured background surface brightness in the CCO 
vicinity. The S/N shows a flat maximum in the range 
1-0" < r < 1.5", with the peak value S/N = 76.5 at 
r = 1.23". To minimize the spectral distortion caused by 
the increase of the PSF size with growing energy (i.e., the 
spectrum extracted from the PSF core is somewhat softer 
than the true spectrum) , we chose the source aperture ra- 
dius of 1.476" (three ACIS pixels), near the larger radius 
of the flat maximum. The MARX simulation shows that 
this circle contains 95.4% of total point source counts and 
94.8% of total energy flux. Because the SNR background 
is not uniform (see Fig. 1), we have to choose the back- 
ground extraction region sufficiently close to the CCO 
but not too close, in order to prevent the background 
"contamination" by the wings of the CCO's PSF. After 
several trials, we chose the 2.46" <r < 3.94" (5-8 ACIS 
pixels) annulus as the background region^. 

We used the CIAO psextract script to extract the 
spectra. The 1.476" radius aperture contains 7016 counts 
in the 0.6-6.0 keV range, of which about 607 counts be- 
long to the background. This corresponds to the source 
count rate of 0.104 ± 0.02 counts s~^ in the source ex- 
traction region (or 0.109±0.02 counts s~^ after correcting 
for the finite size of the extraction aperture). The source 
count rate of ~ 0.033 counts per frame corresponds to the 
pile-up fraction (i.e., the ratio of the number of frames 
with two or more events to the number of frames with 
one or more events) of about 1.6%, a factor of 10-15 

^ To examine the effect of the source and background aperture 
sizes on the results of spectral fitting, we also fit the spectra ex- 
tracted from the source apertures of 1.0" and 2.0" radii, and used 
a few other annular background extraction regions. The spectral 
parameters obtained from those fits are generally consistent with 
those reported below (the differences are within the 2cr uncertain- 
ties). 



lower than in the previous ACIS observations with the 
frame time of 3.24 s. 

We grouped the extracted spectrum with a minimum 
of 30 counts per spectral bin and fit it with several spec- 
tral models. The results are shown in Figures 4 and 
5, and Table 1. The fit with the absorbed power-law 
(PL) model yields the photon index P « 5, much larger 
than P = 1-2 for the X-ray spectra of active pulsars, 
and larger than P = 3-4 obtained from the previous 
observations in which the CCO spectrum was consid- 
erably piled up. The observed flux, F « 6.7 x 10~^^ 
erg cm~^ s~^ in the 0.6-6 keV band (7.0 x 10~^^ erg 
cm~^ s~^ after correcting for the encircled energy frac- 
tion), is slightly lower than that found in the previous 
observations. The flt is marginally acceptable in terms 
of minimum ixt = 1-21 for i/ = 125 dof), but it 
shows excess data counts at lower energies (< 1.2 keV) 
and a deflcit of data counts at higher energies (> 4.5 
keV). This suggests that the true spectral shape is dif- 
ferent from the PL. In addition, the PL fit gives a rather 
large value for the hydrogen column density, A'h,22 ~ 2.8, 
which not only strongly exceeds the total Galactic HI 
column density in that direction (A^hi.22 = .4-0.7; see 
iDickev fc Lockmanll990HKalberla et ai.l2005D but is also 
substantially greater than typical values A'h,22 ~ 1.2-1.3 
(-^11,22 — 1-5 for a maximum value) in the CCO vicinity 
found from the fits of the SNR spectra with the absorbed 
two-compon ent plasma emissi on model in the 1 Ms ex- 
posure data (|Yang et al.ll2008D . 

The fact that the observed spectrum is softer at higher 
energies than the best-flt PL model (even with the large 
P) suggests that a thermal model could provide a better 
description. The flt with the absorbed blackbody (BB) 
model gives^ temperatures kT°° — 0.39-0.41 keV, hydro- 
gen column densities A'h 22 = 1.2-1.4, radii of equivalent 
emitting sphere R°° « (6.8-1.0)i:»3.4 km {R°°/D « 0.27 
km/kpc), and bolometric luminosities LJ^j ~ 3x 10"^^ erg 
s^^. The temperatures are lower, and the radii are larger 
than those obtained from the piled-up spectra from the 
previous observations. The quality of the BB flt is con- 
siderably lower than that of the PL flt. Not only the 
BB flt gives a larger xt (1-54 vs. 1.21), but it shows a 
strong excess of data counts at E > A keV. Therefore, we 
conclude that the BB model gives a poorer description 
of the CCO spectrum than the PL model. 

As the observed spectrum is harder than a BB but 
softer than a PL at higher energies, similar to the spec- 
tra of hydrogen (or helium) NS atmospher e (NSA)^ mod- 
els (jPavlov et al.lll995tlZavhn et al.l |l99d). we fit it with 
a number of such models. An example of such a fit, 
for a hydrogen atmosphere model with a relatively low 

* We add the superscript °° to the temperature and radius to 
show that these parameter are quoted as observed by a distant 
observer, without any corrections for the gravitational effects. The 
temperature, radius and bolometric luminosity "at infinity" are 
connected with those as observed at the NS surface as follows: 
T°° = Tgr, = R/gr, and L^^^ = L^.ig^., where = (1 - 
2GM / Rx?)^^^ is the gravitational redshift parameter. 

" The hydrogen NSA models are dubbed nsa and nsagrav in 
XSPEC. The nsa models are for fixed NS mass and radius, M = 
1.4^^0 and i? = 10 km = 13.06 km, gravitational acceleration 

at the NS surface g = 2.43 X 10^^ cm s~'^); in these models the 
distance is the fitting parameter. The nsagrav models include a set 
of NSA models on a mass-radius grid; therefore, they allow one to 
fit the mass and radius for a given distance. 
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Fig. 4. — Fits of the CCO spectrum extracted from the 1.476" aperture with absorbed one-component spectral models (see Table 1 for 
the fitting parameters). The wabsxNSA panel also shows the background contribution. 



magnetic field {B < G) and NS mass M = IAMq 
and radius i? = 10 km {gr = 0.766), is shown in Fig- 
ure 4. In terms of xt^ the NSA fits are considerably 
better than the BB fit. Most of the NSA fits show an 



excess of data counts at high energies (e.g., 18.3 ± 7.7 
counts in the 4.5-6 keV for the wabsxnsa fit shown in 
Figure 4), but it is not so strong as in the BB fits. Such 
an excess might be caused by the (small) pileup effect 
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(which hardens the count spectrum, as mentioned in 
Sjl]), but our estimates show that only a few such events 
are expected in our observation. The NSA models with 
high magnetic fields {B > a. few 10^^ G) give worse fits 
than the low- field NSA models (in particular, a stronger 
excess at higher energies) because the spe ctra of high- 
field NSA models are closer to BB spectra (jPavlov et alj 
[19951) . In comparison with the BB model, the low-field 
NSA models give a factor of 2 lower effective tempera- 
tures, kT^ = 0.17-0.20 keV, and a factor of ~ 7 larger 
radius-to-distance ratio, R°° /D = 1.3-1.7 km/kpc, for 
M = \AMq and i? = 10 km. For this mass and ra- 
dius, the best-fit distances, D = 7.5-10.1 kpc, substan- 
tially exceed the measured distance to Gas A, suggesting 
a smaller size of the emitting region^°, ^ 4-6 km, 
considerably larger than for the BB fit but still too small 
for a NS radius [and comparable to the Schwarzschild ra- 
dius, Rs = 2.953(M/Mo) km]. This suggests that either 
only a fraction of the NS surface is heated up to X-ray 
temperatures (in which case we should expect pulsations 
in the X-ray emission) or the GGO is not a NS at all (e.g., 
a strange quark star [SQS], whose radius can be consid- 
erably smaller than that of a NS). To explore the latter 
possibility, we fit the GGO spectrum with the nsagrav 
models, which allow one to vary M and R. An example 
of such a fit, in which the mass was fixed at M = O.25M0, 
the distance fixed a,t D = 3.4 kpc, and the radius was 
varied, is presented in Figure 4, for the best-fit param- 
eters R = 5.077 km, T^s = 198 eV, and iVH,22 = 1-57. 
The fit residuals are virtually indiscernible from those 
of the nsa fit (including the slight data excess at higher 
energies), and the (but not Tcfi) and iVH,22 values 
are also almost the same. The corresponding bolometric 
luminosity is Lj^i « 4.2 x 10^^ erg s~\ of which 91% is 
emitted in the 0.6-6 keV band. This fit confirms that 
the observed CCO emission might be interpreted as ra- 
diation emergent from a light-element atmosphere of an 
object more compact than a NS. We will discuss such 
interpretations in more detail in §4. 

If the GGO has a strong magnetic field and its magne- 
tosphcrc contains a large number of energetic electrons, 
these electrons can comptonize the NS thermal emis- 
sion and generate a high-energy tail in the GGO spec- 
trum. This effect can be crudely accounted for by the 
"resonance Gom pton scatte ring" (RCS) X SPEG model 
(|Lvutikov fc Ga vriil 2006; Rea et all [20081 ). The fit of 
the GGO spectrum with the RCS model is statistically 
good {xl = 1-11), but it yields very broad ranges for 
the temperature of the seed BB emission, kT « 0.1-0.2 
keV, and the normalization, w 3 x 10~^-10~^, which are 
strongly anti-correlated (i.e., lower temperatures corre- 
spond to larger normalization parameters). The fit gives 
the resonant scattering optical depth t^^s ~ 1 (at the 
lower boundary of the interval Tjos = 1-10 included in 
the XSPEG set of the RGS models), and the characteris- 
tic velocity (in units of c) of the magnetospheric electrons 
Pt = 0.24-0.35. We see that the seed BB temperature 
is considerably lower than that found from the BB fit, 
which hints that the size of the emitting region is larger 
than BB radius, but, unfortunately, it remains unclear 
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Fig. 5. — Fits of the CCO spectrum extracted from the 1.476" 
aperture with absorbed two-component models (see Table 1 for 
the fitting parameters). The component contributions arc shown 
by dashed lines. 



^" Although the spectral shape depends on the assumed M and 
R, in contrast with the BB spectrum, the dependence is weak, so 
that the radius can be rescaled as ij oc D for crude estimates. 
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how the size could be estimated from the model normal- 
ization (N. Rea, priv. comm.). We should note, however, 
that, since the RCS model itself is substantially over- 
simplified (e.g., it assumes a one-dimensional motion of 
non-relativistic electrons with a rectangular velocity dis- 
tribution along the radial direction), the inferred fitting 
parameters may be quite different from the actual phys- 
ical parameters. Moreover, there is no even circumstan- 
tial evidence that the Cas A CCO has a magnetosphere 
densely populated by energetic particles (such as the non- 
thermal X-ray spectra in rotation-powered pulsars or the 
hard tails seen in magnetars even in quiescence) . There- 
fore, the good fit may simply mean that the RCS model 
just mimics the observed CCO spectrum whose true ori- 
gin may be quite different. 

The quality of the spectral fits can be improved if we 
use various two-component models. Examples of such 
fits are shown in Figure 5 and Table f. For instance, 
the PL+BB model (which is commonly used for the de- 
scription of the X-ray spectra of rotation-powered pul- 
sars, AXPs and SGRs) provides a good fit {x1 = 1.12; 
no excesses or deficits at lower and higher ends of the 
spectrum), but the slope of the PL component is still 
very steep (F « 4.6), the size of the BB-cmitting region 
is still much smaller than the NS radius (« 1.2 km at 
D = 3.4 kpc), and the hydrogen column density is too 
large (7Vh,22 « 2.3). 

The BB-I-BB model (which can be considered as a 
crude model for thermal emission from the surface with 
a non- uniform temperature) gives a good fit (xt = ^■^^)^ 
with fcT,- « 0.3 keV, RZ, « 2 km, L^,,,^,, « 4.1 x 10=^^ 

erg s~^ for the low-temperature (soft) component, and 
fcT^iSrd « 0.6 keV, « 0.25 km, L-,,,,, « 0.6 x 10^3 

erg s~^ for the high-temperature (hard) component (the 
radii and luminosities are for 13 = 3.4 kpc). Although 
the addition of another BB component has resulted in 
a larger size, this size is still much smaller than the NS 
radius. 

Obviously, one could obtain a larger emitting size using 
a two-component NSA model. The NSA-^NSA fit with 
low-field nsa models gives the temperatures ^^^^ « 
0.14 keV and T^^i.^^.^ ~ 0.4 keV. The range of distances 
for the soft component, D = 1.8 -5.9 kpc, includes the 
range of 3.3-3.7 kpc inferred by iReed et all (|1995D for 



the Cas A SNR. Moreover, if we rescale R°°f, 



13.06 



km to i?^ft « 11 km, still consistent with a NS radius, 
it would correspond to the best-fit distance of 3.4 kpc, 
equal to the most probable distance to Cas A. The simi- 
larly scaled size for the hard component is 0.4 km, for 
D ~ 3.4 kpc. For these effective temperatures and scaled 
radii, the bolometric luminosities are LJ^j ^^f^ « 5.7x lO'^'^ 

erg s~^ and L^oi,ha.Td ~ 0.4 x 10^^ erg s~^. Thus, we can- 
not rule out that the CCO is a NS with a nonuniform 
surface temperature, although it is hard to explain the 
nonuniformity without assuming a strong magnetic field 
of a special topology (see §4). 

In all the fits described above, we see some struc- 
ture in the residuals, in the range of ~ 1-3 keV, which 
hints at the presence of spectral lines or photoabsorp- 
tion edges. If real, such spectral features could be ei- 
ther intrinsic features in the CCO spectrum (e.g., similar 
to the absorption lines detected in the spectrum of the 




Fig. 



0.5 1 
Frequency (Hz) 

Power spectrum Z'f{f) for the Cas A CCO. 
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CCO IE 1207.4-5209: ISanwal eTaI][200l ) or they might 
be caused by photoabsorption between the source and 
the observer if the element abundances are different from 
those adopted in the wabs photoabsorption model used 
in the above fits. To explore the latter possibility, we 
applied the vphabs photoabsorption model with variable 
abundances together with the NSA model. Fitting the 
abundances of the elements with photoionization edges 
in the 0.6-3 keV range (Fe, Ne, Mg, Si, S), we found 
that the Fe, Ne, S and Mg abundances were poorly con- 
strained and consistent with the solar values, while the 
Si abundance showed some excess"'^^. Fixing the abun- 
dances of all the elements but Si at their solar values, we 
obtained the Si abundance of 3.2 ± 1.3, without substan- 
tial changes of the other fitting parameters (see Table 
1). The fit shows a marginal improvment with respect 
to the wabsxnsa fit {x1 = 1-19 vs. 1.24) and somewhat 
smoother residuals in the 1.5-2.5 keV range, but it does 
not affect the structure at 1.1-1.4 keV and the apparent 
absorption feature around 2.8 keV. 

We have checked that such features could not be caused 
by inaccurate subtraction of the SNR background, which 
does not show spectral lines at these energies (see the 
panel wabsxnsa in Figure 4). However, our simulations 
of ACIS count rate spectra for continuum spectral mod- 
els show that similar "features" often appear because of 
statistical fluctuations. We also examined the spectra ob- 
tained in the previous, much longer observations of the 
CCO and found no features at these energies for most 
of them^^. Therefore, we conclude that, most likely, our 
observation does not show real spectral features in the 
CCO spectrum. 

3.3. Timing 

To search for the CCO period, we used the arrival times 
for the 6756 events extracted from the 1.23" radius aper- 
ture (which provides maximum S/N; see Figure 3) in the 

Possible additional absorption at the Si-K photoionization 
edge, 1.84 keV, has been noticed by Stage et al. ( 2004) in previ- 
ously observed spectra of the CCO. These authors interpret it as 
caused by material in the Cas A SNR. 

We should mention, however, that the stronger parallel CTI 
effect in those observations might wash out faint spectral features. 
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0.6-6.0 keV band (the estimated background contribu- 
tion is 421 events). With the frame time tfoame = 0.34104 
s, we can search for periods P > 2ttramo = 0.68208 s (fre- 
quencies f < 1.466 Hz). We ca lculated the Zf (Rayleigh) 
statistic ( Buccheri et al.l[r983[ ) as a function of frequency 
(i.e., the power spectrum) in the range 30 ^Hz < / < 
1.466 Hz with the step A/ = l/^Hz « O.OTT^pi^, where 
Tspan = 70, 181 s is the time span of our observation. The 
power spectrum is shown in Figure [6l The highest peaks 
in the power spectrum, — 28.66 and 28.38, are found 
at the frequencies / = 0.814487 and 0.692145 Hz, re- 
spectively (the frequency uncertainty is about 2 /iHz for 
each of the peaks). If, for instance, the highest peak were 
due to the actual pulsations, it would correspond to the 
pulsed fraction p = Pohs{Ns + ^s)^s ^ ~ ['iZl{Ns + 
Nb)Y'^N-^ ± [2{Ns + NB)Y'^Ng^ = 0.098 ± 0.018 
(for a nearly sinusoidal signal), where Ns = 6333 and 
Nb —4:21 are the numbers of source and background 
counts in the source aperture, and Pobs is the pulsed frac- 
tion uncorrected for the background contribution. How- 
ever, the probability of obtaining such a peak by chance 
in a noise spectrum is /maxTspan exp(— Z^/2) = 0.061, 
which corresponds to the confidence level of only 1.87cr. 
Given the low confidence level and the presence of a few 
other peaks of similar heights (see Figure [6|), we have to 
conclude that this observation does not show pulsations 
with a period P > 0.68 s. Following the approach by 
iGroth I (Il975h . the upper limits on the pulsed fraction p 
can be estimated as 13%, 14% or 16% at the confidence 
levels of 95%, 99% or 99.9%, respectively. 

From their anal ysis of the XMM-Newto n observation 
of the Gas A CCO. lMereghetti et al.l(|2002( ) report some- 
what lower values of the pulsed fraction upper limits 
(e.g., p < 13% in the 0.4-2 Hz band at the 3ct level). 
That observation, however, suffered from the strong, 
nonuniform SNR background, whose contribution into 
the source aperture could hardly be measured with a 
sufficient accuracy. Indeed, the fact that the GGO's flux 
reported by Mereghetti et al. is a factor of 3 higher than 
the value we measured with Chandra (see Table 1) sug- 
gests that the background contribution, Nb/INs + Nb), 
is 92% rather than 75% estimated in that paper, and the 
background-corrected upper limit on the pulsed fraction 
is also a factor of 3 higher than reported by Mereghetti 
et al. 

4. DISCUSSION 

Our observation of the Gas A GGO with the Chandra 
AGIS detector in subarray mode has allowed us to obtain 
the most accurate X-ray spectrum, virtually undistorted 
by pileup, to look for extended emission in the high reso- 
lution image of this source, and to search for periodicity 
with periods greater than 0.68 s. The results of this ob- 
servation can be used not only to measure the source 
properties with high precision, but also to put new con- 
straints on its nature. 

The very h igh X-ray-to-opt ical flux ratio (e.g., 
Fx/Fh > 10''; iFesen et al1l2006t ) proves unequivocally 
that the Gas A GGO is indeed a compact object as- 
sociated with the Gas A SNR, not a usual star or a 
background AGN. It is generally believed that a com- 
pact source born in a Type II SN explosion can be either 
a NS or a black hole. If the Gas A GGO were a black 
hole, its X-ray emission could be only due to accretion 



of the ambient matter, in which case we would expect a 
substantial variability (as accretion is an inherently un- 
stable process) and strong emission lines generated in the 
hot accreting matter. However, the GGO's flux has not 
shown stastically significant change s in numerous Chan- 
dra observations (jTeter et al.ll200J ). its X-ray spectrum 
does not show strong emission lines (§3.2), and its X- 
ray luminosity is too high to be explained by accretion 
of SNR material on to the apparently fast-moving GGO 
(jPavlovet al.ll2000f) . Therefore, we conclude that the 
black hole interpretation can be ruled out. 

Most of the known NSs are rotation-powered pulsars 
(~ 1800 such pulsars have been detected in the ra- 
dio, ~ 80 in X-rays, ~ 30 in 7-rays, and ~ 10 in the 
optical). X-ray observations of rotation-powered pul- 
sars usually show hard PL components in their spectra 
(with photon index F ~ 1-2), interpreted as magneto- 
spheric emission, with sharp pulsations. Moreover, vir- 
tually all young, powerful pulsars are accompanied by 
bright PWNe generated by the winds of ultrarelativis- 
tic particles ejected from th e pulsar magnetosphere (e.g., 
iKargaltsev fc Pavlov! l2008l ). As the X-ray spectrum of 
the radio- and 7-ray-quiet GGO is very soft (almost cer- 
tainly thermal; see §3.2), and, more importantly, there 
are no traces of a PWN around the GGO (§3.1), we have 
to conclude that it is not a rotation-powered pulsar. 

The usual (rotation-powered) pulsar activity could be 
quenched in a young NS if its rotation is too slow or 
its magnetic field is too weak to generate strong elec- 
tric fields at and above the NS surface, needed for nor- 
mal pulsar operation. Very high magnetic fields (> 
a few 10^'^ G) can also inhibit the pulsar activity (per- 
haps, by suppressing electron-positron pair production 
beca use of photon splitting i n a superstrong magnetic 
field; 'Baring & Hardmgl [1991) , as we know from obser- 
vations of magnetars. Thus, the Gas A GGO could be- 
long to one of the two known types of young NSs with- 
out pulsar activity — magnetars, which are character- 
ized by superstrong magnetic fields {B ^ 10^^-10^^ G) 
and slow rotat ion (P = 2 -12 s), and "anti-magnetars" 
(Gotthclf & H alpernll2008h . which are NSs with low mag- 
netic fields {B < 10^^ G) and moderately fast rotation, 
whose periods (P a few 0.1 s) have changed very little 
during their lifetime. Also, we cannot exclude the oppor- 
tunity that the GGO represents a new type of a compact 
object, perhaps with an exotic compo sition, such as a 
strange quark star (SOS: Iw"ittenlll984 ). We will discuss 
these options below, with account for our observational 
results. 

4.1. A magnetar? 

A strong support for the magnetar interpretation could 
be provided by detection of pulsations with a period of 
a few seconds, typical for magnetars. We found no evi- 
dence for pulsations, with an upper limit on pulsed frac- 
tion of about 10%-15%, for P > 0.68 s. The nonde- 
tection, however, does not prove that the GGO is not 
a magnetar bec ause some of them have a n even lower 
pulsed fraction ([Woods fc Thomps"oiill2006l ). 

The shape of the GGO's X-ray spectrum resembles 
those of magnetars in the Chandra energy band. For 
instance, the PL+BB model satisfactorily describes the 
spectra of both the G GO and magnetars. Howeve r, ac- 
cording to Table 4.1 in lWoods fc ThompsonI ([2006f) . the 
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magnetar BB temperatures (0.4-0.7 keV) inferred from 
such fits are somewhat higher than the temperature of 
the BB component in the CCO's spectrum, 0.30-0.38 
keV, while the magnetar photon indices (F = 2-4) are 
somewhat smaher than the photon index of the CCO's 
PL component, F = 3.6-5.2. The fitting parameters 
obtained from the RCS fit s are also close to those ob- 
tained by iRea et all (|2008f ) for some AXPs. However, 
the CCO's X-ray luminosity, Lx ^ 4 x 10'^^ erg s~^, is 
lower than those of magnetars, Lx ~ lO'^^-lO'^^ erg s~^. 
It would be very interesting to compare the spectra of 
the CCO and magnetars at i? > 10 keV, where magne- 
tars show very hard PL spectra (F = 1-1.5 for the total 
[pulsed + nonpulsed] emission), but it would require hard 
X-ray instruments with high spatial resolution to sepa- 
rate the CCO's emission from the much brighter Cas A 
SNR's emission. 

If the CCO were a magnetar, we would expect bursts 
and flares. However, we have not seen any significant 
flux variations in the 10 years of Chandra observations. 
In fact, if in the last ^ 40 years the CCO had produced a 
flare on a scale we see in SGRs, it would likely has been 
detected by an all-sky mon itor in hard X-rays or soft 7- 
rays. There was a claim bv lKrause et al.l ()2005[ ) that the 
apparent motions of Cas A filaments with tangential ve- 
locities close to the speed of light, observed with Spitzer 
at 24 /im, could be interpreted as an infrared echo of a 
strong flare from the CCO ~ 60 years ago, but th is in- 
terpretation has been r etracted bv lKim et al.l (|2008l ) and 
iDwek fc Are^ (poM ). 

One might speculate that the CCO does possess a 
superstrong magnetic field, perhaps hidden in the NS 
interior, but, being so young, it has not yet devel- 
oped the properties we see in "mature" magnetars. 
Such an interpretation is in line with the hypothesis by 
[Bhattacharya & Soni (2007), who suggest that the su- 
perstrong magnetic field in a very young magnetar is 
originally concentrated in the NS core, and it will require 
> 10'^ yr for this field to reach the crust and give rise to 
the magnetar activity. It is hardly possible to confirm 
or reject this interpretation of the CCO by direct obser- 
vations (unless a burst or a flare is detected). In the 
framework of this hypothesis, the X-ray emission from 
the CCO should be thermal, with luminosity and tem- 
perature perhaps somewhat higher than those of "usual" 
NSs due to additional heat that might be released in 
the NS interiors by dissipation of the "hidden" super- 
strong magnetic field, and its spectrum should be best 
described by an NSA model. The observed spectrum gen- 
erally does not contradict this hypothesis, at least if we 
believe that the NS surface temperature is nonuniform, 
but this would not be a unique interpretation. The inter- 
pr etation will perhaps lo o k mor e plausible if the scenario 
bv lBhattacharva fc Sonil (|2007f ) is supported by detailed 
studies of established magnetars (i.e., AXPs and SGRs). 
To conclude, the magnetar interpretation of the Cas A 
CCO does not seem very likely, but it cannot be ruled 
out, especially the assumption that the CCO is an "im- 
mature" magnetar. 

4.2. An anti-magnetar? 

In addition to the Cas A, CCOs in several other SNRs 
are known, with sim ilar properties (jPavlov et al.l [2003 , 
I2004t Ide Luca1l2008f) . Those CCOs also show thermal- 



like X-ray spectra, albeit with temperatures and lumi- 
nosities lower than those of the Cas A CCO, and they 
also do not show pulsar or magnetar activity. It seems 
natural to explore the possibility that the nature of the 
Cas A CCO is the same as of the other CCOs, just the 
Cas A CCO is younger. 

Unlike the Cas A CCO, pulsations have been dis- 
covered in three CCOs, with periods in the range 
of . 105-0.424 s (IZavliri et al.l [2000I : iGotthelf et all 
I2OOI IGotthelf fc Halpernl \2Qm . Surprisingly, the up- 
per limit s on period derivatives turned out to be 
very low (IGotthelf fc Halpernll2007l 120091 : iHalnern et all 
|2007| ). corresponding to spin-down luminosities, E = 
4:Tr'^IPP~^ , lower than the X-ray luminosities (which 
means that the X-ray emission is not powered by NS 
rotation). To explain the fact that the corresponding 
lower limits on spin-down ages are much larger than the 
SNR ages, Tgd = P/2P » Tsnr, one has to assume 
that the initial spin periods of these pulsars were very 
close to their current values, Pq = P{\ — Tsnr/tsci)^^^- 
As the upper limits on the magnetic fields^ , B = 
3.2 X 10i9(PP)i/2 < a few X 10" G, turned out to be 
much lower than those of magnetars (and even of young 
rotation-powered pul sars), these objects were d ubbed 
"anti-magnetars" by IGotthelf fc Halpernl (|2008D . who 
also suggested that the Cas A CCO may be a member of 
this class of NSs. 

The current data, including our new observation, do 
not contradict the assumption than the Cas A CCO is 
indeed an anti-magnetar. The nondetection of period is 
not conclusive because the shortest period we would be 
able to detect is longer than the periods of the known 
anti-magnetars, and the upper limit on pulsed fraction is 
not low enough to exclude weak pulsation even at these 
long periods (e.g., the pulsed fractions of the CCOs in 
the PKS 1209-51/52 and PuppisA SNRs are about 9% 
and 11%, respectively). The Cas A CCO spectrum is 
well described by the same spectral models as those of 
anti-magnetars, with similar spectral parameters. Simi- 
lar to the anti-magnetars, the spectra are thermal (or at 
least they contain a strong thermal component), but a 
single-component BB model does not give good fits. A 
single-component NSA model gives acceptable fits, but 
the radii of equivalent emitting sphere are smaller than a 
reasonable NS radius. Two-component models, such as 
BB-I-BB, give very good fits, but the NS surface layers 
at so high temperatures and low magnetic fields are ex- 
pected to be gaseous (plasma) atmospheres, whose spec- 
tra are different from BB because photons with different 
energies emerge from layers with different temperatures. 
Therefore, a more realistic description of the thermal 
emission from anti-magnetars should be based on NSA 
models rather than the BB model. 

As the magnetic field of an anti-magnetar is low, and 
the spin period is relatively long, accretion of SN debris 

13 For the best-investigated CCO, IE 1207.4-5209 in the 
PKS 1209-51/52 SNR, the magnetic field, Bcycl = 6 X I0^°gr^ 
G w as estimated from the absorption features IjSanwal et al.l 
120021) interpreted as h armonics of the electron cyclotron frequency 
dBignanu "etani2003h . A very close value of the magnetic field, 
_Ba(j = 9.7 X 10^" G, has been recently estimated from the period 
deriva tive determined with the aid of new XMM-JVewton observa- 
tions llPavlov et al.ll2009bl) . 
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soon after the SN explosion is possible. If there is even a 
very small fi'action of hydrogen or helium in the accreted 
material, the emission emerges from a hydrogen or he- 
lium NS atmosphere because of the gravitational strati- 
fication. If the field is much lower than 10^^ G, it makes 
no effect on the NSA X-ray emission, so that the use of 
low-field NSA models is justified. If the field is around 
10" G (as in IE 1207-5209), cyclotron absorption fea- 
tures are expected in the Chandra energy band. This 
would require the use of NSA models with the cyclotron 
energy in the soft X-ray range (currently unavailable in 
XSPEC), but the continua of such models are rathe r close 
to the low-field NS spectra (jSuleimanov at al.ll2009( ). The 
relatively low fit quality and too small size of the emit- 
ting region in the single-component NSA fit suggests 
that the NS surface temperature is not uniform. As 
we do not know the temperature distribution over the 
surface, we have to use the simplest approximation of 
two regions with different temperatures and areas (e.g., 
a hot spot and a colder bulk surface). Our fit with the 
NSA-I-NSA model has given a very good fit, with tem- 
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11 km and i?hard ~ 0.4 km, for D = 3.4 
kpc. The bolometric luminosities of the two components 
are L^,, 5.7 x 10^3 erg s'^ and L?^,^ ~ 0.4 x 10^3 
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s^", assuming isotropic emission. Since the bolometric 
luminosity and the temperature of the soft component 
are consistent w ith the predic tions of various NS cool- 
ing models (e.g.. lTsurut"alll998l ) for the Gas A GGO age, 
the soft component can be interpreted as emission of the 
thermal energy stored in the NS interiors. However, it is 
not easy to explain the origin of a small hot spot at the 
NS surface. At such low magnetic fields, the anisotropy 
of thermal conductivity, which may lead to a nonuni- 
form temperature distribution (e.g., the magnetic poles 
are hotter than the equator because the conductivity is 
higher along the magnetic field), is negligible, and, more- 
over, this effect could not produce such a small hot spot 
even in a superstrong magnetic fi e ld for the dipole ge - 
ometry (" Perez- Azon'n et all [20061) . iPavlov et all ()2000r ) 
speculated that the temperature anisotropy could be due 
to a nonuniform chemical composition over the NS sur- 
face (e.g., there might be hydrogen "islands", perhaps 
formed by accretion onto the magnetic poles soon af- 
ter the SN explosion, on the iron surface that is colder 
than the islands because the thermal conductivity of the 
degenerate NS envelope is proportional to the inverse 
nucleus charge, Z~^). However, this interpretation looks 
somewhat artificial, and it has not been supported by ob- 
servational data (e.g., heavy-element spectral lines in the 
GGO spectrum). Alternatively, the hot spot(s) could be 
heated by ongoing slow accretion from a residual debris 
disk onto the magnetic pole(s), but there are no obser- 
vational indications, such as flux variability or emission 
spectral lines, of such accretion. Another way to explain 
the origin of small hot spot(s) is to assume the pres- 
ence of a very stron g toroidal magnetic field, i? > 10^^ 
G, in the NS crust (jPerez-Azorfn et al.l [20061) . Such a 
field could screen the heat flux from the NS interior to- 
ward the surface in a broad range of magnetic latitudes 
around the equator, leading to a large temperature con- 
trast between t he poles and the rest of the NS surface (see 
Figs. 8 and 9 in lPerez-Azorm et al.l[200l l. Although the 
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Fig. 7. — Mass-radius region compatible with the distance to 
the ceo, obtained from the wabsx nsagrav fits. The black strip 
corresponds to D = 3.3-3.5 kpc, while the white strips leftward and 
rightward of the black strip correspond to D = 3.1-3.3 and 3.5- 
3.7 kpc, respectively. The dotted lines correspond to i? = Rg and 
-R = 1.5-Rg, where Rg is the gravitational (Schwarzschild) radius; 
the solid lines, 314 = cms~'^ = 10 and 0.1, correspond to 

the maximum and minimum gravitational accelerations for which 
the nsagrav models are available. The M(R) curves for two models 
of strange quark matter equation of state, SQMl and SQM3, are 
shown by dashed curves. The dash-dotted curves are the loci of 
constant R°° (note that R.°° R. for M 0). 

toroidal field would not penetrate into the NS surface 
layers and would not affect the NS spin-down (which is 
determined by a much lower poloidal [e.g., dipole] compo- 
nent in this model), the hypothesis about the presence 
of a very strong toroidal component in anti-magnetars 
may look somewhat artificial. In a sense, within the 
framework of this hypothesis, anti-magnetars become 
similar to the "immature m agnetars" in the scenario by 
iBhattacharva fc Sonil (|2007D . 

To conclude, the interpretation of the Gas A GGO as 
a young anti-magnetar seems very plausible, but there 
is no an unequivocal explanation for the nonuniformity 
of the surface temperature. This problem, however, is 
pertinent to the established anti-magnetars, too. The 
anti-magnetar interpretation of the Gas A GGO could 
be confirmed by a measurement of the GGO period with 
a small period derivative, which would require dedicated 
observations^"* with a timing resolution considerably bet- 
ter than 0.1 s. A firm proof of a low surface magnetic 
field could be provided by detection of two or more spec- 
tral features associated with cyclotron harmonics, which 
would require deep observations with an observational 
setup similar to that used in our observation. To re- 
solve the problem of nonuniform temperature of anti- 
magnetars and understand the temperature and field 
distributions over the NS surface, phase-resolved spec- 
troscopy, interpreted with the aid of NSA models with 
appropriate magnetic fields, would be most helpful. 

4.3. A strange quark star? 

As we have seen in §3.2, reasonably good fits to the Gas 
A GGO spectrum are provided by the NSA models with 

A 300 ks observation with the Chandra HRC-S detector, 
aimed at the search for the CCO period, has been carried out re- 
cently (2009 March; PI D. Chakrabarty), but no results have been 
published yet. 
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a uniform effective temperature, but these fits imply a 
small emitting area, much smaller than the visible area of 
a NS. To infer the allowed range of masses and radii of the 
ceo, assuming that it is covered by a low-field hydrogen 
atmosphere, we have fit the observed spectrum with the 
nsagrav models on a grid in the R-M plane. Each of 
such fits yields the distance to the source, through the 
model normalization. This allows one to plot the lines 
of constant distance (and strips corresponding to ranges 
of distances) in the R-M plane. Figure 7 shows such 
strips around the line corresponding to D — 3.4 kpc , 
the most probable distance to Cas A (jReed et al.lll995[ ): 
the black strip {D = 3.3-3.5 kpc) is enveloped by the 
strips for D = 3.1- 3.3 and 3.5-3.7 kpc. We see that 
the maximum allowed radius R of the compact object 
is about 5.5 km, while the apparent radius R°° < 6.5 
km, and the object's mass cannot exceed O.SMq. This 
R-M domain is inconsistent with any NS equation of 
state, but it is consistent with low-mass SQS models. 
The M{R) curves for two equation s of state of strange 
quark matter, SQMl and SQM3 (Lattimer fc PrakashI 
1200 It) , shown in Fi gure 7, cross the R-M domain allowed 
by the nsagrav fits at i? « 4.5-5.5 km; they correspond 
to the allowed mass ranges M « 0.2-0.3 and 0.1-0.2 Mq, 
respectively. The effective temperature and bolometric 
luminosity, « 2.1 MK and Lj^j « 4 x 10^3 erg s^^, 
do not change significantly within the strips. 

Thus, the emission from the CCO could, in prin- 
ciple, be interpreted as emerging from a SQS atmo- 
sphere which could form above a normal-matter crust 
([Glcndcnning & Wcbcr 1992). The main problem with 
this interpretation is the very low star's mass. For in- 
stance, if the SQS is formed as a result of a phase tran- 
sition in a NS, we would not expect the SQS mass to 
be lower than w IMq. One might speculate that the 
SQS was formed directly in the SN explosion, but we are 
not aware of SN models that would result in such a low- 
mass compact remnant. On the other hand, we cannot 
be sure that such a scenario is impossible. Therefore, 
we believe that the SQS interpretation, however exotic 
it looks, should not be dismissed until rejected by fur- 
ther observations. A possible way to check it would be 
very deep observations in the near-IR (more suitable than 
the optical because of the high extinction) , which could 
detect the Rayleigh- Jeans tail of thermal emission from 
the bulk of the NS surface. If the small region associated 
with the CCO X-ray emission is, in fact, just a heated 
area on the NS surface, while the rest of the surface is 
substantially colder, the near-IR observation could de- 
tect this colder emission and prove that its size is much 
larger than that seen in X-rays. The current deepest 
near-IR limit, infer red from HST NICMOS observations, 
is H > 24.6 (.Fesen et al.|[2006i) . but the emission from 
the NS surface is expected to be even fainter, i.e., deeper 
observations are required to detect it. The SQS interpre- 
tation would also be called into question if future X-ray 



observations show this object to be an anti-magnetar. 
(Although anti-magnetars might, in principle, be SQSs, 
the more conventional NS interpretaion is preferable un- 
til proven otherwise, according to the Occam's razor.) 

4.4. Summary 

To summarize, we believe that the anti-magnetar in- 
terpretation of the Cas A CCO, similar to the three other 
CCOs proven to be anti-magnetars, is currently the most 
plausible. In this interpretation, the CCO's emission 
with the bolometric luminosity Lj^j ^ 6 x lO'^^ erg s~^ 
emerges from a NS atmosphere with nonuniform effec- 
tive temperature, relatively low magnetic field, B < 10^^ 
G, and very slow spin-down. The origin of the tempera- 
ture nonuniformity is not clear, as well as for the other 
anti-magnetars; it might be caused by a much stronger 
toroidal magnetic field in the NS crust. The most di- 
rect way to check the anti-magnetar interpretation is to 
measure the CCO's period and period derivative and/or 
detect spectral features that could be interpreted as har- 
monics of the electron cyclotron frequency. 

There is no observational evidence of the Cas A CCO 
being a magnetar similar to the currently known SGRs 
and AXPs. We, however, cannot rule out the possibil- 
ity that it is an immature magnetar, whose ultrastrong 
magnetic field is hidden in the NS interiors and does not 
make a strong effect on the observable properties. If this 
is the case, the Cas A CCO may eventually turn into 
an "ordinary" magnetar, but we consider this hypothesis 
rather speculative. 

Finally, we cannot firmly reject the possibility that the 
Cas A CCO is a SQS, with a radius of !=s 5 km and 
a mass ^ O.SMq, covered by a normal-matter crust and 
an atmosphere with kT^ « 0.18 keV and Lj^j w 4 x 10^^ 
er s~^. It remains unclear how the SQS with so low mass 
could be formed. Future observations will help check this 
interpretation. 
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Fitting parameters for different spectral models 
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Note. — The fits arc for the source spectrum extracted from the 1.476" radius aperture, background taken from the 2.26"— 3.94" 
annulus. The errors of fitting parameters are given at the 90% confidence level for one interesting parameter. 

^ Hydrogen column density in units of lO^^ cm~2.t> photon index for PL spectra."^ Effective temperature as seen by a distant observer, 
in units of eV. PL normalization = photon spectral flux at 1 keV, in units of 10"^ photons cm"^ keV^l.'^ Radius (as seen by a 
distant observer) in units of km, and distance in kpc. In the nsa model, the NS radius and mass are fixed at i? = 10 km and M = 1.4Mq 
{R°° = 13.06 km), while the distance is the fitting parameter. ' Observed (absorbed) energy flux in the 0.6-6 keV band, in units of 
10~13 ergs cm~2 s~^, calculated with the cflux model in XSPEC.^ Unabsorbed energy flux in the 0.6-6 keV band, in units of lO^^^ ergs 
cm~2 s , calculated with the cflux model in XSPEC. The corresponding luminosity is Lx = 1-38 x s~i, for D = 3.4 kpc. 

^ NSA model for fixed M = IAMq and R = 10 km.' The fit is for fixed M = 0.25Mo. The distance (normalization) was fixed, while 
the radius i? is a fitting parameter {R = 5.08lQ'gg km).j The result is for variable Si abundance; its fitted value is 3.2lj^'j.'' Fit with the 
resonance Compton scattering model, for fixed Tres = 1-0. In addition to A^h a-nd kT, the fitting parameters are Pt = 0.3310:02 and norm 
= 0.71J00 X 10-4. 
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